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EVIDENCE OF AN INITIALLY MAGNETICALLY DOMINATED OUTFLOW IN GRB 080916C 
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ABSTRACT 

The composition of gamma-ray burst (GRB) ejecta is still a mystery. The standard model invokes 
an initially hot "fireball" composed of baryonic matter. Here we analyze the broad band spectra of 
GRB 080916C detected by the Fermi satellite. The featureless Band-spectrum of all five epochs as 
well as the detections of > 10 GeV photons in this burst place a strong constraint on the prompt 
emission radius i?^, which is typically > 10^^ cm, independent on the details of the emission process. 
The lack of detection of a thermal component as predicted by the baryonic models strongly suggests 
that a significant fraction of the outflow energy is initially not in the "fireball" form, but is likely in 
a Poynting flux entrained with the baryonic matter. The ratio between the Poynting and baryonic 
fluxes is at least ~ (15 — 20) at the photosphere radius, if the Poynting flux is not directly converted 
to kinetic energy below the photosphere. 

Subject headings: gamma rays:bursts — gamma rays: observations — gamma rays: theory — plasmas — 
radiation mechanisms: non-thermal — radiation mechanism: thermal 



1. INTRODUCTION 

In the classical model of Gamma-ray bursts (GRBs), 
an initially hot "fireball" is composed of photons, elec- 
tron/positron pair s , and a sinall amount of baryons 
JPaczyhskiJ 119861 iGoodman I Il986t IShemi fc Pir"an1 
[19901. This fireball is soon accelerated to a relativis- 
tic speed under its own thermal pressure. Due to 
the existence of baryons, a significant fraction of en- 
ergy is converted into the kinetic energy of the ejecta 
(jMfezaro s et al. 1 119931 : [Piran et al. Ill993l ). The rest of 
energy is still stored in the form of photons, which 
escape the system as the fir eball becomes transpar- 
ent at the photosphe re radius (|Meszaros fc ReesllBool : 
iMeszaros et al. Il2002f ). The ejecta then coast with a rel- 
ativistic speed without significant radiation until reach- 
ing a larger radius when a fraction of kinetic energy 
is dissipated into heat an d radiation in internal shocks 
ees fc Meszaros I Il994l ). The internal shock model 
has the advantage of interpreting GRB variabili ties, but 
suffers the low radiation effi ciency problem ([Kumar I 
[19991 : [Panaitescu et al. I [1999[ ). Recently, an analysis 
of the prompt emission data of several GRBs suggests 
that the internal shock model is disfavored by the dat a 
([Kumar fc McMahon I [2008L [Kumar fc Naravan I [2001 ). 
An alternative model invokes a dynamically important 
magnetic field. Within such a Poynting flux domi- 
nated model, the observed GRB emission is powered 
by dissipation of the magneti c energy within the ejecta 



([Usov |[1994l: [Thompson 
[Meszaros fc Rees H1997F 



T99llV lahakis & Konigl 2003; 



Lvutikov fc BlandfordM2003,) . 



Until recently, it has been difficult to diagnose the GRB 
composition from observational data. The recent detec- 
tion of the broad-band featureless Band-function spectra 
as well as the very high gamma-ray emissio n (> 10 GeV) 
from GRB 080916C by the Fermi satellite ([Abdo al.\ 
[2009[ ) provides a unique opportunity to diagnose the 
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GRB composition^. Below we will use the opacity argu- 
ment (§2) and the photosphere argument (§3) to argue 
that the ejecta of GRB 080916C must contain a signifi- 
cant fraction of energy that is initially not in the "fire- 
ball" form, but in a Poynting flux entrained with bary- 
onic matter. 



2. MODEL INDEPENDENT EMISSION RADIUS 
CONSTRAINT 

Both the Large Area Telescope (LAT) and the Gamma- 
ray Burst Monitor (GBM) on-board Fermi have detected 
GRB 080916C. The time-dependent spectral analysis re- 
veals a series of featur eless smoothly-join t broken power 
law ( "Band-function" , [Band et al. 1[1993[ ) spectrum cate- 
gorized by a peak energy Ep and two asymptotic ph oton 
spectral power law indices a and /3 ([Abdo et al. [ |2009). 
These spectra cover 5-6 orders of magnitude in energy, 
from - 10 keV to ~ (1 - 10) GeV. The highest pho- 
ton energy reaches 13.2 GeV (in the t ime interval "d" 
defined in Fig.l of 'Ab do et aU ([2009D 'l. At a redshift 
z = 4.35 ± 0.15 (Grein er et an[2009[ ). this burst is the 
most energetic GRB known to date, with an isotropic 
gamma-ray energy ~ 9 x lO^'* erg. 

The gamma-ray spectrum is expected to have a pair 
cutoff feature at large energies due to the compactness 
argument, i.e. the optical depth for two photon pair pro- 
duction (77 — > e"'"e~) may reach unity above a critical 
cutoff energy Ecut- Within the internal shock model, the 
pair cutoff energy (or the lack of it) , together with the ob- 
served variability time scale, can be used to constr ain the 
bulk Lorentz factor P of the outflow (iB aring fc Ha rding 1 
|l997l : lLithwick fc Sari |[200ll : [ Abdo etal. 112009 ). In some 
models, the observed variability time scale may not re- 
flect that of the central engine ([Lvutikov fc Blandford I 
[20031: [Naravan fc Kumar [[20091 ). So more generally the 

An 18- GeV pho ton was detected earlier from GRB 941017 
ijHurlev et al. I [19941 1. However, it was significantly delayed, not 
associated with the prompt emission. The bandpass of prompt 
emission observation was not wide enough to constrain the thermal 
emission component. Hence, it could not be applied to directly 
constrain the GRB emission radius and composition. 
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pair cutoff energy can be expressed as a function of two 
independent parameters, the bulk Lorentz factor T and 
the g amma-ray emission radi us Rj from the central en- 
gine ijGupta fc Zhang 1 12008[ ). For each cutoff energy 
Ecut, one can define a threshold energy Etd above which 
the photons with this energy can interact with the pho- 
tons at Ecut to produce pairs. This threshold condition is 
defined by {Etd/1 MeV){E,ut/l MeV)> 0.25[r/(l + z)]^. 
The expression of the pair production optical depth T-y-y 
depends on the relative location between Etd and the 
Band-function break energy. For Ep{a — f3)/{2 + a) < 
Etd < Ecut (which is generally satisfied for this analysis), 
the optical depth of a photon with the observed energy E 
can be coasted into a simple form (derived from Eqs.(13) 
and (14) of Gupta & Zhang 2008, or from Eqs. (3-4) of 
Lithwick & Sari 2001 with 6t absorbed into the expres- 
sion of and with the cosmological correction factor 
(1 -|- z) properly taken into account) 

T^^{E) = 



1 



r 



-1-/3 V^eCV -R7 VI + ^ 

(1) 

where dz = {c/Hq) dx/y^Q\ + ilm(l + a;)^ is the co- 
moving distance of the GRB, z is the redshift of the GRB, 
and TOe, c and or are the fundamental constants elec- 
tron mass, speed of light, and Thomson cross section, 
respectively. Here the energies, i.e. E and meC^, are 
both in units of "keV". The parameter /g (in units of 
ph • cm^^ • (kcV)^^^'^) is such defined that the observed 
photon fluence spectrum above the break is written as 
M{E) — fpE^- I n terms of the fitting parameters of 
lAbdo et aL \ (l2009h . /o can be expressed as 



fo = A-AT 



Ep{a-P) 



(2 + a) 



a-/3 



exp(^-a)(100keV)" 



(2) 

where AT is the (observed) time interval taken to per- 
form the Band-function fit, and A is the Band-function 



normalization in units of ph • cm 



•keV" 



Fi- 



nally, the coefficient C{f3) in Eq.(l) is a function of /3, and 
we adopt the approx imation C {{3) ~ (7/6)(— /3)~^/^/(l — 
(3) (jSvensson II1987D to perform the calculations^. 

The lack of a spectral cutoff feature suggests 
r-y^(£'niax) < 1, where i?max is the maximum energy of 
the observed photons. Using this condition, one can de- 
rive the r — R^ constra i nts fo r the 5 time intervals (a-e) 
defined bv lAbdo et aTl (|2009D ( using the values in their 
Table 1, with AT defined by the time ranges listed in 
the table). By doing so, we have implicitly assumed that 
the emission location and Lorentz factor for a particular 
time interval AT essentially remain the same. Figure 1 
shows the five critical lines (solid) in the F — i?^ space. 
The allowed parameter regimes are above the lines. For 
all the 5 time intervals the allowed emission radii are 
all large. This model-independent conclusion regarding 
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Fig. 1.— The T-R diagram of GRB 080916C. The constraints on 
Fty are displayed in color thick lines, above which are the allowed 
parameter spaces for each corresponding epoch. The photosphere 
radii Rpi^ are displayed as color dashed lines, with the same color 
convention. The internal shock radii i?ig = V'^cSt°^ /(I + z) for 
St°^ = 0.5, 2 s are the black dashed lines. 

Swift observations reveal a commonl y seen steep-decay 
phase in the early X-ray afterglow (|Tagliaferri et aL\ 
l2005f ). which suggests that the GRB prompt emission 
region is separated from the emiss^ion region of the after- 
glow (the external shock) (|Zhang al. Il2006[ ). Within 
the "internal" models, the expected dissipation radius 
spans in a wide range, from the photosphere radius (typ- 
ically at Rph ~ (10^^ — 10^^) cm) to slightly smaller 
than the deceleration radius (~ 10^^ cm). The de- 
rived large emission radius R^ is inconsistent with the 
photosphere model (see §3 for detail), but is consistent 
with that expected from the m agnetic dissipation model 
(|Lvutikov fc BlandfordI l2003f ). The "internal shock" 
model invokes an emission radius i?is ^ T^c5t (where 
5t is the variability time scale of the central engine). 
The T — Ry constraint of this model is plotted in figure 
[1] as black dashed line s for two values of the variability 
time scales discussed in lAbdo et al. I (|2009( ). 5t = 0.5, 2 s. 
The internal shock site is allowed by the data. The con- 
strained internal shock radii are at least Ris ^ 10^^ — 10^^ 
cm for the five epochs. The required Lorentz factors are 
at least 500-1000 (see Fig.l). However, there are further 
constrains on this model, as discussed below. 

3. CONSTRAINTS FROM THE PHOTOSPHERIC 
(THERMAL) EMISSION 

Besides the energy dissipation region (internal shock or 
magnetic dissipation site), the fireball "photosphere", at 
which the fireball becomes transparent during the expan- 
sion, i s another e mission s ite of GRBs (Me szaros fc Reesi 
2000t iRees fc Meszarosi l2005t [Thompson al. I l2007l: 



Ghisellini et al. 1 120071 : iLazzati et al. I I2009D . The pho- 
tosphere radius Rph is defined by the electron scatter- 
a large R^ is consistent with the results de rived from ing optical depth r^^ = n'arl^ = 1, where n' and A' 



other independent methods for o t her GRBs (iLvutikov I 
2OOI iKumar fc McMahmTI [2008t [Racusin et al. 1 120081: 



Kumar fc Naravan |[2009l: IShen fc Zhang |[2009l) . 



^ IGupta fc Zhang I 1I2OOI) takes C(/3) = 3/8(1 - /3), while 
ULithwick fc Sari 1120011) ta kes C(/3) = 11/180. For a typical value 
/3 = - 2, the C(/3) value of ISvenssoiTI (1987) and Gupta fc Zhang 



(2008) are consistent with each other, but that of .Lithwick fc Sari I 
1 2001 ) is smaller by a factor of ~ 2. 



are the electron number density and width of the ejecta 
shell in the rest frame comoving with the ejecta. Be- 
low we will assume a pure baryonic flow (as expected 
in the fireball model) to derive the expected photo- 
sphere spectra. For a continuous baryonic wind (which 
is suitable to describe GRB 08091 6C given the smooth 
hghtcurves reported in Fig.l of lAbdo et oTI (HqOOi)) 
with a total wind luminosity L^,, the photosphere ra- 
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dius can be written as (|Meszaros fc Reesl[2000t 

er I 

120081) Rph ^ {L^aTRo/STrmpC^r])^/^ for Rph < Rc, and 
Rph ~ Ly^UT / ^T^TTipC'V^ri for Rph > Rc- Here, nip is 
the proton mass, rj — L^/Mc^ is the dimensionless en- 
tropy of the baryonic flow, and Rc ^ Rq ■ min(?7,7y*) is 
the radius at which the ejecta wind reaches the "coast- 
ing" phase. The critical dimensionless entropy is ?]* = 
{LwUt /itTTmpC^ RqY^^ , and Rq = cSt is the radius where 
the ejecta is emitted from the central engine. The coast- 
ing Lorentz factor is F = ry for Rph > Rc, and T = rj^ 
for Rph < Rc- The photosphere radii for the 5 epochs of 
GRB 080916C are presented in Fig.l by the color dashed 
lines. 

It has been suggested that the observed Ep of GRBs 
may be the therrnal peak of the photosphere em ission 
(jRees fc Meszaros I I2005t [Thompson et al. I I2007D . and 
that the Comptonization of the thermal e mission may 
lead t o a n on-thermal spectrum above Ep (jPe'er et al. I 
[2005l[20nl . However, for GRB 080916C, Figure 1 sug- 
gests that the observed gamma-ray emission is not from 
the photosphere, since i?-y ^ Rph- Under the high- 
compactne ss condition, a second pair photosp here may 
be form ed (M eszaros et al. ll2002t[Kobavashi et al. 1120021: 
iPe'er fc Waxman 1 120041) 7 but its radius is also much 
smaller than R^ inferred from the data. Another sugges- 
tion is that the observed GRB spectra are the superposi- 
tion of a thermal component (photosphere) that defines 
E^ and a non-t hermal power law component ( Ryde l[2005t 
iRvde fc Pe'er ir200&) . This model is based on the data in 
a relatively narrow spectral range of BATSE, and when 
extrapolated to high energies, would violate the observed 
featureless Band-function significantly. The theoretically 
expected thermal peak energy « 50 keV (see below) is 
we ll below the o bserv ed Ep in the 5 epochs discussed 
bv lAbdo et~ar\ (|2009( l. Clearly this model by itself is 
insufficient to interpret the data of this burst. 

The initial wind luminosity of the fireball is at least 
the observed gamma-ray luminosity (assuming radi- 
ation efficiency 100%), i.e. L^, > L^. Such an outflow, 
during the expansion, must have released a residual ther- 
mal emission at the transparent (photosphere) radius. 
The luminosity of this thermal componen t is high, and 
can be written as ([Meszaros fc R.ees l[2000h 

_ J itu, ?7 > 7?*, Rph < Rc, ,o\ 

"^"^ - \ L^{r,/Tj,f/^ rj < ij,, Rph > Rc. 

The expected temperature of the black body compo- 
nent emerging from the photosphere is thus T°jl = 
(L^/47ri?2ca) 1/4(1 + z)-^ for Rph < Rc, and = 

{L^/^TTRlcaf'\Rph/Rc)-^'''{l + z)-^ for Rph > Rc, 
respectively, where a is Stefan-Boltzmann's constant. 

Assuming a blackbody spectrum (other modified spec- 
tral shapes would not change the conclusion), we plot 
in Fig. 2 the lower limit of the expected photosphere 
spectrum for the internal shock model in the baryon- 
dominated outflow (by taking = Lj) and compare 
it with the observations. For clarity, only the prediction 
for the epoch (b) are plotted (red thick-dashed curve). 
This plot corresponds to rj > r/^ = 850, hence Lth = L^, 
Rph 0^ Rc and = T°l^^^^ = 50 keV. It is appar- 
ent that such a component strongly violates the observa- 
tional data, which is a featureless Band-function. 



One can further argue that the baryonic model in gen- 
eral does not work. In Fig. 2, we present the predicted 
lower limits of the photosphere spectrum for two addi- 
tional (more conservative) temperatures T^^ = 10, 1 keV. 
These correspond to either a coasting Lorentz factor 
rj = 470, 200 (for the same 6t°^), or a central engine vari- 
ability time scale much longer than what is observed (for 
the same rj), which may be relevant to a fireball emerg- 
ing from a stellar envelope in the coUapsar model. These 
models give a much larger Rph, and hence, a much lower 
T°jl. Since St"'' = 0.5 s has been observed (|Abdo et al. I 
l2009f ). the models that invoke a longer central engine 
variability time scale require the unconventional assump- 
tion that the obser yed variability time scale is not that of 
the central engine (|Naravan fc Kumar II2009I ). Nonethe- 
less, the predicted thermal spectra of these cases (red, 
thin-dashed curves) also strongly violate the observa- 
tional constraints. 

This analysis strongly suggests that the initial wind 
luminosity was not stored in the "fireball" form at 
the base of the central engine. Since there is no 
other known source of energy that can be entrained in 
the ejecta, we identify the missing luminosity as the 
Poynting flux luminosity, which is not observable before 
strong magnetic dissipation o ccurs at a much larger ra- 
dius. It has been s uggested (jZhang fc MeszarosTl2002l : 
iDaigne &: Mochkovi tch 2002|) that the photosphere ther- 
mal component can be much dimmer if the outflow is 
Poynting flux dominated. In order to hide the bright 
thermal component, one can pose a lower limit on the 
ratio between the Poynt ing flux and the ba r yonic flux, 
a = Lp/Lb,. Following I Zhang fc Meszaros I (|2002D . we 
define L„ = Lb + Lp = (1 + a)Lb- Assuming no dissi- 
pation of the Poynting flux below Rph, the above pho- 
tosphere derivations can be modifled by replacing L^j 
by Luj/{1 + a). One can then derive the required min- 
imum a value that can "hide" the photosphere thermal 
component. For the specific internal shock model re- 
quired by the variability time scale (the red thick-dashed 
curve), one requires at least cr ~ 20 at the photosphere to 
make the photosphere emission unobservable (red thick- 
dotted curve in Fig. 2). At such a high-cr shocks are very 
weak and c annot power the bright gamma-ray emission 
(I Zhang ^Fl Cobayashi 2005), so unless a can reduce sig- 
nificantly from the photosphere to the internal shock ra- 
dius (e.g. Spruit et al. 2001) but without generating a 
dominant emission component in the observed spectrum, 
the internal shock model is not a viable mechanism to 
interpret this burst^. For the other two unconventional 
cases, similar minumim values of cr = 15 — 20 at photo- 
sphere are required in order to obtain consistency with 
the observed spectra. We can therefore conclude that 
the ejecta of GRB 080916C cannot be a pure "fireball", 
but must store a significant fraction of energy initially in 
a Poynting fiux. 

4. SUMMARY AND DISCUSSION 

In this letter, we show that the observed feature-less 
Band-function broad band spectra of GRB 080916C ob- 

Another possibility would be that the Poynting flux is con- 
verted to kinetic energy directly without mag netic dissipation dur- 
ing the acceleration phase (e.g. IVlahakis fc Konigl n2003i V The 
photosphere brightness of this model is not studied in detail. 
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Fig. 2. — T he observed Band-function spectra for the five epochs 
IIAbdo et al. 2009) (color solid) and the predicted lower limits of 
the photosphere spectra (red dashed) for different parameters for 
the epoch (b) within the framework of the baryonic fireball models. 
Red, thick-dashed curve: the internal shock model with 5t°'' = 0.5 
s, corresponding to T°j^ = 50 keV; red, thin-dashed curves: for 
= 10, 1 keV. The suppressed photosphere spectra are plotted 
by red, dotted curves, with the required cr values marked. 

served by the Fermi satelhte (|Abdo et al. I l2009t ) pose 
two interesting constraints on the GRB prompt emis- 
sion mechanism. First, the detection of high energy, 
> 10 GeV, photons places a strong model-independent 
r — constraint on the GRB ejecta (Fig.l), which pre- 
cludes the photosphere {Rph <^ R-y) as the emission site. 
Second, the non-detection of a bright thermal compo- 
nent expected from the baryonic models puts a strong 
constraint on the composition of the fireball: the flow 
should be dominated by a Poynting flux component with 
a minimum a « 15 — 20 at the photosphere in order 
to account for the observed spectra (Fig. 2), as long as 
the Poynting flux energy is not directly converted to the 
kinetic energy of the flow below the photosphere. 

Within such a picture, the observed bright gamma-ray 
emission may be powered by dissipation of the Poynt- 
ing flux energy within the outflow. A Poynting flux 
dominated flow favors synchrotron radiation as the emis- 
sion mechanism of the observed gamma-rays. This 
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is consistent with independent modeling of the burst 
(jWang et al. Il2009h . It also suggests a weak synchrotron 
self-Compton component, which is not observed from the 
data. This model requires that the GRB central engine 
launches a collimated magnetic bubble from the central 
engine (black hole or millisecond magnetar), which prop- 
agates through the st ellar envelope withou t reducing the 
a value significantly (| Wheeler et al. |[2000D . After escap- 
ing the star, the ejecta is quickly accelerat ed under its 
own magnetic pressure ('Mi zuno" et al. 11200 9^. At the ra- 
dius where photons become transparent (photosphere), 
the bulk of the wind energy is not in the thermal form. 
The cr-value may reduce above photosphere. However, in 
order to have internal shocks, the magnetic energy needs 
to be "quietly" dissipated below the internal shock radius 
without significant emission. More naturally, the mag- 
netic field may be globally dissipated in a cataclysmic 
manner to power the observed emission. This model 
does not demand a high a at the deceleration radius, 
since a of the flow would decrease significantly in the 
prompt emission region due to intense magnetic dissipa- 
tion. In any case, the a value at the deceleration ra- 
dius can be of order unity or higher, which gives a weak 
or moderately bright reverse shock emission com ponent 
(jZhang fc Kobavashi |[2005HMimica et al. 11200^ . 

This analysis applies to GRB 080916C. Similar anal- 
yses can be carried out for more bursts co-detected by 
Fermi LAT/GBM in the future to determine whether 
other GRB outflows are also magnetically dominated. 
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